Dipole matrix element approach vs. Peierls approximation for optical conductivity 
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We develop a computational approach for calculating the optical conductivity in the augmented 
plane wave basis set of Wien2K and apply it for thoroughly comparing the full dipole matrix element 
calculation and the Peierls approximation. The results for SrVOs and V2O3 show that the Peierls 
approximation, which is commonly used in model calculations, works well for optical transitions 
between the d orbitals. In a typical transition metal oxide, these transitions are solely responsible 
for the optical conductivity at low frequencies. The Peierls approximation does not work, on the 
other hand, for optical transitions between p- and rf-orbitals which usually became important at 
frequencies of a few eVs. 

PACS numbers: 71.27.+a,71.1.Fd 



Much of our knowledge about solid state systems 
comes from their response to small electro-magnetic per- 
turbations. A broad range of techniques has been devel- 
oped to probe the nature of ground states in elastic scat- 
tering experiments and the excitations in inelastic scat- 
tering or absorption experiments. It is usually a com- 
bination of several experimental techniques as well as 
theoretical calculations which allow us to draw a com- 
plete picture of a given material. Among those, the opti- 
cal spectroscopy plays an important rolei, complement- 
ing the photoemission spectroscopy (PES) which is eas- 
ier to access and interpret in most theories. Probing the 
particle-hole excitations averaged over the Brillouin zone, 
the optical spectra contain a different and less detailed 
information about the system than angle-resolved pho- 
toemission spectra. The main asset of the optical spec- 
troscopy, however, is its robustness: Unlike PES, it does 
not suffer from surface effects. Moreover, unlike trans- 
port measurements, the optical conductivity is not crit- 
ically affected by impurities or disorder: Optical transi- 
tions cannot simply disappear, but can only be shifted to 
different energies, which is expressed by the sum-rule for 
optical conductivity^. 

Calculations of optical spectra from first principles are 
well established within the effective non-interacting elec- 
tron theories^ for weakly correlated materials such as 
the local density approximation (LDA)^ to the density 
functional theory. The many-body perturbation theory 
on the GW level^ and its two-particle extensions us- 
ing the Bethe-Salpeter equatior*^ have been successful 
in describing the excitonic physics in semiconductors. 
The situation is different in the field of strongly corre- 
lated electron systems. Although the optical measure- 
ment on these materials proved very useful for investi- 
gation of metal-insulator transitions or mass renormal- 
ization effects, material specific theoretical investigations 
are rather rare»l"— . This is perhaps not surprising given 
the fact that calculation of one-particle spectra is already 
a formidable challenge. 

In the past decade the dynamical mean-field theory 



(DMFT)i^"— combined with first-principles bandstruc- 
tures (LDA-I-DMFT) ^^'^^ showed considerable power to 
describe correlated materials. This theory allows for an 
accurate description of the local (intra-atomic) dynam- 
ics, while the inter-atomic effects are treated on the static 
mean-field level. Importantly, DMFT is not restricted 
to a particular energy scale and thus allows for the si- 
multaneous description of quasiparticles on meV scale 
and atomic excitations on the eV scale, which is crucial 
to capture the spectral weight transfers in the optical 
spectra. First optical calculations with DMFT were per- 
formed for the single-band Hubbard model. It was shown 
that the local approximation of DMFT leads to vanishing 
of the vertex corrections to the optical conductivityi^"— . 
This means that the electron and hole created in the 
process of optical excitation behave independently and 
thus the Green's function of the electron-hole pair is a 
product of two one-particle Green's functions. This is 
not necessarily true in multi-band models, except for the 
case of degenerate bands. However, the dipole selection 
rules at optical frequencies typically forbid creation of an 
electron-hole pair on the same atom and thus the vertex 
corrections may be neglected also in this case, an approx- 
imation we also adopt throughout this work. Note that 
for inelastic x-ray scattering experiments 'optical' tran- 
sitions with finite momentum transfer allow formation of 
strongly bound local electron-hole pairs, excitons. The 
vertex corrections in this case are substantial. A typical 
example are the crystal-field d-d excitations deep in the 
optical gap of transition metal oxides; see Refi^ for a 
V2O3 calculation. 

While the formal framework for calculating the opti- 
cal conductivity within the above approximations is well 
established, the numerical implementation poses several 
challenges: i) k-space integration, ii) determination of the 
optical transition amplitudes and inclusion of states in a 
broad energy window, iii) evaluation of optical spectra 
for real frequencies, which is an additional problem aris- 
ing for particular numerical techniques, such as quantum 
Monte-Carlo simulations, used to solve the DMFT equa- 



tions. Different strategies for dealing with tfiese issues 
are possible. In this article, we present an implementa- 
tion based on the Wannier functions formalism and a 
direct calculation of the transition amplitudes from the 
one-particle wave functions. We compare our results to 
the so-called Peierls approximation, which relies on the 
k-dcrivatives of the effective low-energy Hamiltonian of 
the systems considered and discuss their relationship. 
We analyze specifically two well-known correlated oxides, 
SrVOa and V2O3, as archetypes, and compare our results 
to the available experimental data. 

The outline of the paper is as follows: In Section HI we 
give details on the LDA-I-DMFT calculation. In Section 
III Al the technical details of the dipole matrix element 
caculation are discussed. Section fll Bl discussed the rela- 
tionship to the Peierls approximation, which is popular 
for lattice models. Sections IIII Al and llll BI present the re- 
sults for SrVOs and V2O3, respectively. Finally, Section 
IIVI summarizes the main findings. 



I. LDA+DMFT WITH WANNIER ORBITALS 

The DMFT equations are naturally formulated in 
terms of fermionic creation and annihilation operators 
on a lattice, a formulation which assumes an underlying 
set of localized orthogonal orbitals. Our starting point 
are the LDA Bloch states V'nk and corresponding band- 
energies Snk calculated with the full-potential linear aug- 
mented plane waves (LAPW) program Wien2k<^i^ De- 
pending on the specific material considered, we choose 
an energy window defined by the lower and upper band 
indices nmin l£ n < rimax and transform the states ipnk 
from this window to real-space Wannier orbitals^^ local- 
ized around lattice sites R: 
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where C/(k) are unitary matrices defined throughout the 
Brillouin zone and A^k is the number of k-points. Us- 
ing wien2wannier— and wannier902^ the matrices C/(k) 
that lead to maximally localized Wannier functions are 
found. Construction of the single-particle part of the ef- 
fective Hamiltonian is completed by rotation of the LDA 
Hamiltonian into the Wannier basis 
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Finally the on-site interaction is added to the Hamil- 
tonian. Another important input data required for the 
DMFT calculation, arc the local Coulomb repulsion pa- 
rameters which define the term to be added to H^^^,(k) 
in the Wannier basis: the intra-orbital local repulsion 
[/, inter-orbit local interaction V and the exchange pa- 
rameter J. In principle, this input should be computed 
from the underlying LDA data, with constrained LDA^ii 
or constrained random phase approximation^. However, 



since identifying U, V, and J for SrVOa or V2O3 is 
not the aim of this work, we adapted these values from 
the literature^""— iS^. In the case V2O3 we have chosen 
a slightly lower value of U than in Refsi^""— , which 
ensures the best agreement with XAS^^ and optical 
experiments^, according to considerations reported in 
RefjH. It is important to notice though that the Coulomb 
parameter as well as the DMFT itself in general depend 
on the chosen basis set of Wannier function, which espe- 
cially becomes important if the choice of Wannier orbitals 
is not as straightforward as in the case of SrVOa or V2O3 
below. In both cases the actual DMFT calculation has 
been done in the t2g subspace. For the optical conductiv- 
ity, the DMFT t2g Green function was then supplemented 
with the LDA Green function for the other orbitals. 

Once the effective LDA Hamiltonian is set up, the 
DMFT equations are solved numerically using quantum 
Monte Carlo simulations^ with an imaginary time dis- 
cretization of At := 0.1 eV"^ for SrVOa and At = 0.125 
eV~^ for V2O3, respectively, to obtain the one-particle 
self-energy which serves as the many-body input for eval- 
uation of the optical conductivity. Similar LDA+DMFT 
approaches based on augemted plane waves can be found 
in Ref.^^>^. 

From the LDA-I-DMFT self energy T,{iujm) for a tem- 
perature T at the Matsubara frequencies Wm = TT{2m + 
l)r, we obtain S on real frequencies via the procedure 
described in Ref42i: Starting from the imaginary-time 
Green's function G{t) measured in Monte-Carlo, the k- 
integrated spectrum A{uj) is calculated by the maximum 
entropy method (MEM, see Ref4S). Afterwards, the local 
Green's function Gmem['^) for real frequencies is found 
by applying Kramers-Kronig relations. Finally, we fit 
S(w) such that the Green's function obtained by direct k- 
summation, i.e. G{uj) = 1/A^k X)k['*^~^W'(k) — S(w)]~"^, 
matches the one from the maximum entropy method. 



II. LINEAR RESPONSE FOR THE OPTICAL 
CONDUCTIVITY 

A. Dipole matrix element approach 

The regular part of the optical conductivity is obtained 
via the standard Kubo's formula in linear response^ 

a^piu) = hm Re (^ [ dt e'^' ([ja(q, t), j;3(-q, 0)]) 

(3) 

where V is the unit cell volume, ia{^.^t) is the q- 
momcntum paramagnetic current in a direction, ft = 1 
and q -^^ is the dipole approximation. Expressing Eq. 
([3]) via the Lindhardt bubble and Matsubara formal- 
ism (thus omitting vertex corrections, consistently with 
the discussion above) , we include internal degrees of free- 
dom describing optical excitations^: initial (final) fre- 
quency w,(aj -|- J7), reciprocal vector k and N = Umax — 



nmin orbital degrees of freedom participating in optical 
transitions. Altogether, we obtain the following expres- 
sion for the real part of the optical conductivity 
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where / is the Fermi function, wi^(k) = t/(k)w"(k)[/+(k) 
are rotated matrix elements of the momentum operator 
-y^„(k) = -i(-0nk|VQ|V'mk)/me, 1 <,n,TO < A^, the el- 
ementary charge e = 1, and A(k, w) = — ImG'(k, cij)/7r is 
the generalized spectral function with the Greens func- 
tion 



G(k,cj) = [{uj + ii)l- Hwi^i) - S(w)]" 



(5) 



Here, // denotes the chemical potential, H\\r(}^) G C^^^ 
the (non-interacting) Hamiltonian in the Wannier or- 
bital basis and S(aj) G f^^^xJV ^^^^ ggj£ energy from the 
LDA-f DMFT calculation. 

For an efficient and accurate k-quadrature of Eq. (j4]) , we 
use a tetrahedral-mesh integration. To resolve regions in 
k-space with larger integration error we adaptively refine 
the tetrahedra in these domains. Furthermore, the sym- 
metry operations of the unit cell are applied such that 
the integrand of Eq @ has to be evaluated only at k- 
points within the reduced wedge of the Brillouin zone. 
The computation of the momentum matrix elements 
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which are in the following also denoted as dipole ma- 
trix requires their evaluation in term of the underlying 
LAPW basis set^. It is thus (to our knowledge for the 
first time) possible to combine a full potential LAPW 
dipole matrix with Wannier-functions-based DMFT al- 
gorithms for the computation of transport and optical 
properties (for different approaches see e.g. Refs4^i^). 
Note that the surveyed workflow is not limited to the use 
of a DMFT self energy S(a;), but can be easily general- 
ized for other, even k-dependent self energies S(k, w). In 
such cases, however, the inclusion of vertex corrections 
to the bubble term becomes usually necessary^. 
In addition to transitions within Hilbert space of the low 
energy model, the present approach also allows inclusion 
of higher energy bands. This can be achieved by enlarg- 
ing the transformation matrices C/(k) and, consequently, 
the Hamiltonian Hwi}^) 
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with diagonal En (k),i?„ ^(k) ~ e„k for n < rimin {n > 
nmax)- Note that though the U and Hw are block- 
diagonal, the corresponding dipole matrix V(k) = 
Z^(k)w"(k)W+(k) is not. Inserting iY, ?^, V into Eq.® and 
(O, we thus also take transitions between the Wannier 
orbitals and Bloch states outside of the low energy model 
into account. 



B. Peierls approximation 

For many-body calculations of lattice models, it is 
common practice to determine the optical conductivity 
by the Peierls approximation(PA)i^. The PA approxi- 
mates the group velocities directly from the hopping el- 
ements and, for non-Bravais lattices, from the atomic 
positions in the unit cell^iiSi^. If one wants to go be- 
yond the PA, however, one needs to know the underly- 
ing continuum description for calculating the dipole ma- 
trix elements. The idea of the PA is a gauge transfor- 
mation of the electro-magnetic potential A which disre- 
gards the inner orbital structure (an orbital will get a 
different gauge factor at different positions) and assumes 
a single gauge factor which only depends on the lattice 
site. This is refiected in a modified hopping amplitude 
iRm;R'n -^ tRrn;R'n exp[iA(R - R')/c}i^i^ between sites 
R and R'. 

In the following we discuss the corrections to the PA, 
emerging from the exact continuum description in the 
Wannier orbitals basis, cf. Ref. — . Using the operator 
identity — — V = [Hq , r] , where Hq is the one-particle 
part of the Hamiltonian we can write the momentum 
matrix element as 

{Wrnk\iy\w,n'k) = ^ y^e«k(R'-R) {w^^\[Ho, r]\Wm'R') 

m _ /v ^ — ^ 
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The first term equals Vk-H^(fc), the PA, and can be ob- 
tained without explicit knowledge of the orbital, e.g. from 
an empirical tight-binding Hamiltonian. The remaining 
two terms can be further analyzed by noting that the 
Wannier functions form a complete eigenbasis of Hq. 
Hence we have 



(Wmk|iV|u'™'k) = Vki?mm'(k) +Cmm'(k), (10) 
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Let us first discuss the corrections for a single atom in 
the unit cell. These can be classified as follows 
(i) Intra-atomic dipole transitions: Terms in Eqs . ([9|) and 
dni) withR = R' yield together {w„iii\[Ha,r]\wm'ii), i.e., 
the atomic-dipolc elements with the only difference being 
that Hq is the one-particle Hamiltonian of the solid and 
not of the atom. These local transitions generally require 
different angular momenta for m and m' orbitals and are 
hence at a higher energy. They cannot be described by 
the PA which only considers a single gauge factor for the 
atom or site. 

(ii) Dipole transition mediated hopping: For R' = 
L 7^ R, the first term of Eq. ([Tl]) consists of a hop- 
ping integral {wm'R\Ho\wpRi) and a local dipole tran- 
sition (wpR'|r — R'jwni/R,/). This is similar as intra- 
atomic dipole transitions, however now we obtain a k- 
dependence which was absent for (i). Note, the same 
is obtained for the second term of Eq. pi[l in the case 
R = L y^ R'. 

(iii) Inter-atomic dipole transitions: For R = L 7^ R', the 
first term of Eq. pTJ) consists of a local Wannier matrix 
element {'Wmii\Ho\'WpB^) and an inter-atomic dipole tran- 
sition (wpRJr — R'jWm'R,'). A similar term with a minus 
sign is obtained for the second term of Eq. fTT]) in the 
case R' = L 7^ R. If the orbitals are locally orthogonal, 
only the local on-site energies survive, and we only get 
a contribution if there is a crystal field splitting of the 
orbitals. 

(iv) Further corrections arise if all lattice positions R, 
L and R' are different in Eg. pT]) . In this case we have 
a combination of an inter-atomic dipole element and a 
hopping term. 

If the orbitals are more localized, i.e., exponentially de- 
caying between the atoms, both the hopping element and 
the inter-atomic dipole element are affected by this expo- 
nential suppression. Hence the terms (iv), which contain 
two such exponentials, are more strongly suppressed than 
the hopping amplitude itself (which enters (ii) and (iii) 
as well as the PA) since it only contains one exponential 
factor. The leading term in the "localized" limit is (i), 
which only involves local transitions. 

From these general considerations, the PA appears a 
rather unjustified approximation. In fact, even in the 
limit of more localized orbitals only the terms (iv) get 
suppressed. However, in specific cases of interest the PA 
may be justified. For instance, terms (i)-(iii) become only 
relevant if the orbitals are (a) affected by a large crys- 
tal field splitting or (b) of a different angular momen- 
tum, which typically also means large excitation energies. 



Hence for transitions below this energy, e.g. the Drude 
peak, PA is expected to work, at least for sufhciently 
localized orbitals. 

The situation becomes a bit more involved in the 
case of several atoms in the unit cell. Tomczak and 
Biermann^iiS showed that the PA has to be generalized 
to include the hopping terms between the atoms in the 
same unit cell, which are absent in Vkffmm'(k). How- 
ever, also in this case, the same correction terms (i)-(iv) 
as discussed above remain. 



III. RESULTS 
A. SrVOa 

Due to its simple cubic (perovskite) lattice struc- 
ture and 3d^ electronic structure, SrVOa has been em- 
ployed as a testbed for ab initio calculations such as 
LDA-I-DMFT. There is, on average, a single d electron re- 
siding in three degenerate t2g bands that cross the Fermi 
energy Ep. These t2g bands are well separated by a gap 
from the oxygen p bands below and the e^ orbitals above. 
This situation makes the electronic structure of SrVOs 
particularly simple. 

The photoemission spectra^"— show a well devel- 
oped lower Hubbard below Ep band and a pro- 
nounced quasiparticle peak around Ep; an upper Hub- 
bard band is found, on the other side, in x-ray ab- 
sorption experiments!^ The quasiparticle peak is renor- 
malized (narrowed) by a factor of about 2 compared 
to the overall LDA t2g bandwidth^! This is in good 
agreement with LDA-I-DMFT calculations^i2^ in which 
the interaction parameters have been determined from 
constrained LDA calculations^. Essentially the same 
one-particle spectrum has also been obtained in sub- 
sequent LDA-I-DMFT calculations (among others see 
Refs. 27,32,5i- _56^. g^j^j^ various Wannier function projec- 
tion schemes have been tested for this prototypical ma- 
terial (among others see Refsi^"— ). SrVOa is also the 
materials where kinks in strongly correlated electron sys- 
tems, abstain from any external bosonic degrees of free- 
dom and anti-ferromagnetic spin-fluctuations have been 
discovered f^i^ Similar structures can also be identified 
in angular resolved photoemission spectra<^ As the opti- 
cal conductivity averages (integrates) however over differ- 
ent k-points, such fine structures are hardly discernible 
in this physical quantityi^ Experimentally, the optical 
conductivity shows a Drude peak and additional fea- 
tures above 2 eV when transitions between Hubbard and 
quasiparticle peak become relevant (among others see 
Refs,^!^). 



1. Spectral properties 

The LDA density of states (DOS) for SrVOa used in 
our analysis can be found in Fig. [T] (top panel), where 



the three partial DOS contributions V-t2g, V-e^ and O- 
p are highlighted (we use a = 3.84 A as unit cell lat- 
tice parameter for the cubic perovskite). From the LDA 
data, we obtained three different Wannier projections: 
First, just the V-t2g manifold was mapped onto three 
Wannier orbitals (in the following abbreviated as PI). 
Second, we also included the two additional bands with 
predominant V-e^ character and thus describe the full 
V-d manifold (P2) . Finally, we also take into account the 
O-p bands which leads to a basis consisting of 14 Wan- 
nier functions (P3). 

In Fig. [1] middle and lower panel, we plot the 
LDA-I-DMFT spectra computed with the Wannier ba- 
sis sets PI and P2, respectively. The parameters where 
adapted from Ref.— : local intra-orbital Coulomb repul- 
sion U = 5.05 eV, local inter-orbital repulsion V = 3.55 
eV, and local exchange J = 0.75 cV. Compared to the 
LDA DOS, the t2g partial density of states is renor- 
malized and the formation of lower and upper Hubbard 
bands can be observed as correlation effects are taken 
into account within the DMFT framework. In the case of 
P2, where all 5 V Wannier orbitals are included in the 
DMFT, the Cg orbitals remain completely unoccupied as 
in LDA leading to negligible correlation effects in these 
two orbitals (see lower panel of Fig. [1]). We thus restricted 
the LDA+DMFT analysis for lower temperatures to PI, 
where only the t2g orbitals are described within DMFT. 



Comparison of the dipole matrix elements approach and 
the Peierls approximation 
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FIG. 1. (Color online) Non- interacting partial density of states 
of SrVOs (results abbreviated as LDA) compared to DMFT 
spectra at T = IIGOA'. Two DMFT basis sets were employed: 
first, the three orbital t2g basis with parameters {U, J, V) = 
(5.05, 0.75, 3.55) eV; second, the entire V d manifold with the 
same interaction parameters for all orbitals. 



In Fig. [51 our main results for the optical conductiv- 
ity of SrVOs are summarized. We compare four differ- 
ent calculations for the (isotropic) optical conductivity 
a computed via Eq.Q with the experimental data from 
Refi^: The uppermost panel shows a computed by use of 
the LDA Greens function ([5]), where we fixed the broad- 
ening by setting S = —O.OAi [eV] in Eq.®, and em- 
ployed the dipole matrix ^ as group velocities. The sec- 
ond panel of Fig. [21 visualizes the optical conductivity 
a computed with the same Green function G'(k, w), but 
with the Peierls approximation VkHwOf) for the group 
velocities (we are neglecting for this calculation the intra- 
unit-cell contributions as introduced by Tomczakiii) . For 
the lower two panels, we inserted the DMFT self energy 
into the formula for the Greens function ([5|). In partic- 
ular, the third and the fourth panel of Fig. [51 show the 
LDA-I-DMFT results for a calculated with the dipole ma- 
trix and the Peierls approximation as group velocities, 
respectively. Note that the effect of taking a different 
temperature in the experiment (T = 290 K) and in the 
calculations (T = 460 K) is expected to be limited since 
in this temperature range SrVOa does not show a notable 
change in the electronic structure. The main consequence 
of lowering the temperature T = 460 — >■ 290 is the de- 
creased electron-electron scattering within the coherent 
part of the electron spectrum which eventually leads to 



the Drude peak to become more pronounced while the 
inter-band contributions remain essentially unchanged. 
In our analysis of the results, let us start investigating 
the qualitative effect of correlation on the optical spec- 
tra, i.e. comparing the upper two panels of Fig.[2lwith the 
two lower ones. The renormalization of the t2g manifold 
surveyed in Fig. [1] leads to a smaller Drude weight in the 
LDA+DMFT panels of Fig. [21 and to a suppression of the 
prominent peak of a around 3.5 eV predominately stem- 
ming from transitions from the occupied O-p manifold 
to the unoccupied section of the V-i2g orbitals. The sup- 
pression of these two features is also seen in experiment. 
Additionally, the DMFT optical spectra in the lower two 
panels of Fig. [2l show the formation of a small satellite at 
^ 2 eV originating from transitions from the lower Hub- 
bard band of the t2g orbitals to the unoccupied part of 
their coherent spectral peaks. 

Comparing the dipole matrix approach with the Peierls 
approximation, i.e. the first with the second and the third 
with the fourth panel in Fig. [21 indicates both reproduce 
low energy transitions in a similar way. Since the Drude 
peak in this material stems from intra-band excitations of 
the t2g bands, this implies that the Peierls approximation 
is sufficient to describe optical transitions in SrVOs as 
long as only well localized orbitals are participating. The 
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FIG. 2. (Color online) Optical conductivity of SrVOs calcu- 
lated with dipole matrix elements and the Peierls approxima- 
tion, respectively, compared to experimenli^. 



case is different for the 0-p -^^-tig transitions, where a 
deviation in the range between 4 eV and 13 eV is clearly 
visible. Here, the dipole matrix matrix approach appears 
to be superior and a is much closer to the experimental 
results than the Peierls approximation, especially for the 
LDA-(-DMFT optical spectra (see third and fourth panel 
of Fig. [2]). The reason for this behavior can be understood 
taking into account the more non-local nature of the 0-p 
orbitals and the deficiency of the Peierls approximation 
to describe optical transitions therein quantitatively cor- 
rect. 

In addition to deviations compared to the choice of group 
velocities, both LDA-I-DMFT results for the optical spec- 
tra deviate with experiment around 3.5 eV at the onset of 
the 0-p — >'V-i2g transitions. Since LDA seems to describe 
this onset more accurately, we think that the reason for 
this behavior deduces from the fact that, including only 
the tig orbitals in LDA-I-DMFT, we did not consider a 
double counting correction shifting the tig orbitals rela- 
tive to the 0-p orbitals. A more complete approach would 
consider the 0-p within LDA+DMFT on the level of the 
Hartree approximation taking into account the double 
counting corrections more accurately. Then, the change 
of the tig orbitals within LDA+DMFT would eventu- 
ally shift the -p states to lower energies correcting the 
energy distance between the onset of the 0-p manifold 
to the (now renormalizcd) peak in the tig orbitals back 
to the LDA level. 
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FIG. 3. (Color online) Optical conductivity of SrVOs com- 
paring LDA, DMFT and experiment. The LDA and DMFT 
results were computed by the dipole matrix element ap- 
proach (bottom). The top panel shows the sumrule Neff{Qc) 
from Eq. (|12|l for a of the lower panel (the experimental data 
are from Ref.^. 



3. Sumrule analysis 

An important aspect associated to the theoretical and 
experimental study of the optical spectroscopy is the 
analysis of the associated f-sum rule^. This is a direct con- 
sequence of charge conservation, stating that the integral 
over all frequencies of the optical conductivity is always 
proportional to the total electronic density ritot = Ntot/V 
of the system 



lim — 2 / duj a(uj) = lim 



fia-^oo Tre 



Nt, 



V 



V 



(12) 

The importance of the f-sum rule, however, goes well be- 
yond the verification of the charge conservation in LDA 
or LDA-I-DMFT calculations of optical spectra. The val- 
idation of Eq. p^ in theoretical calculations as well as 
in experiment represents a rather academic but delicate 
issue, as it involves very different energy scales (corre- 
sponding to optical transitions involving valence and core 
states). For further details about this issue, we refer the 
reader to Refsi^i^. 
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FIG. 4. (Color online) Optical conductivity of V2O3 (a-phase, 
metallic) at T = 460 K calculated with dipole matrix ele- 
ments [a: LDA, b: DMFT for V-t25 with (U,J,V) = (4,0.7,2.6) 
eV] and the Peierls approximation [c: LDA, d: DMFT for V- 
t2g with (U,J,V)=(4,0.7,2.6) eV], respectively, compared to 
experiment (dashed line taken from Ref.— ). 



More specific information can be extracted by the anal- 
ysis of so-called partial or restricted optical sum-rules. 
They correspond to consider just a portion of the fre- 
quency integral in Eq. 1121 a typical case being a finite up- 
per cut-off ^c^ and how this partial integral changes as a 
function of external parameters (e.g., temperature, mag- 
netic field, etc.). This provides usually very important 
information about the energy balance associated, e.g,. 
with a phase transitions, as it has emerged from many 
experimental^ and theoretical analysis^ of integrated 
optical spectroscopic data of high-temperature supercon- 
ducting cuprates, and most recently, by analyzing2^ the 
non-Slater nature of the antiferromagnetic phase in the 
optical spectra of LaSrMnO^S. 

An example for the application of Eq. (fT2)) is reported 
in the upper panel of Fig. |31 where the growth of A^e// 
with increasing frequency up to Vic = 20 eV is shown 
for the case of SrVOs (there are 19 valence electrons in- 
cluded in our calculation). While a detailed analysis of 
the restricted sum rules for SrVOs goes beyond the scope 
of this work (for the analysis of the restricted sum rule in 
V2O3, see, e.g., RefjS), when comparing the integrated 
LDA and LDA+DMFT optical spectra of Fig. [3l we can 
note, for the latter case, a slight decrease of the values 
of Neff in the low-frequency region, which reflects, ev- 
idently, a correspondent reduction of the electronic mo- 
bility due to electronic correlations. At higher frequency, 
however, the LDA electron density value has been recov- 
ered within an accuracy of about 3 %. 



B. V2O3 

Vanadium scsquioxide V2O3 has been subject of a con- 
siderable interest in condensed matter physics since the 



early Seventies (see e.g. Reffii)^ as it represents one of 
the most evident realization of the Mott-Hubbard metal- 
to-insulator transition (MIT). In fact, V2O3 can be rela- 
tively easily doped with Cr or Ti, and its phase diagram 
displays a clear first order transition between a para- 
magnetic metallic (PM) state (at low concentration of 
Cr, or for Ti doping) and a paramagnetic insulating (PI) 
at a higher level of Cr doping. Such a first order MIT, 
which emerges from a (simultaneous) lower temperature 
structural and magnetic transition and ends up at higher 
Ts with a second order (critical) endpoint, is completely 
isostructural: The high-T paramagnetic phases of Cr^;- 
V2-a;03 are always associated with a corundum crystal 
structure. 

The experimental evidence of the MIT in V2O3 has 
been accumulated, first for static quantities (e.g., the 
dc resistivity) and -ar a later time- for spectral func- 
tions (PESl^, ARPES^, XAS^ili, etc.). In this paper, 
however, we focus on infrared-optical spectroscopy^^^i^ 
only, which is a bulk sensitive technique in comparison 
to photoemission, and -contrary to XAS- includes impor- 
tant information about the itinerant part of the electronic 
properties of strongly correlated electron systems. In op- 
tical spectroscopy measurements at room T, the cross- 
ing of the MIT upon Cr-doping is clearly reflected in 
the abrupt disappearance of the (weak) Drude peak in 
the in-plano2S optical conductivity (t{ijj) with the open- 
ing of a sizable spectral gap. Further important informa- 
tion has been also extracted from the temperature^^ and 
pressure^ dependence oi a{uj): The former has provided 
a clear indication of a strong interplay between small 
changes of the lattice parameters and electronic proper- 
ties, while the latter (together with XAS measurements 
of the V K-pre-edge) has proven the inconsistency of the 
long-standing assumption of equivalence of doping-level 
and applied pressure in the phase diagram of V2O3. Also 
to be mentioned are very recent optical measurements^ 
performed in the most "intriguing" region of the phase- 
diagram, i.e., right across the MIT first order transition 
line: The combined analysis of optical data and photoe- 
mission on a microscopic scale has demonstrated the for- 
mation of insulating islands embedded in the PM phase 
in the metallic side of the MIT. The formation of such is- 
lands, growing in size when the transition is approached, 
can be put -to some exert- in analogy with the nucleation 
processes due to impurities in a standard liquid-gas tran- 
sition: In the case of V2O3 the impurity would be likely 
provided by the lattice distortions^^ due to the Cr- sub- 
stitutions. 

From the theoretical point of view, the problem to 
be analyzed consists of a systems with two electrons in 
the three 3d-t2g (i.e., correlated) bands of the V atom 
at the Fermi level. The t2g basis further splits in one 
aig and two e^^ local orbitals (separated by 0.2-0.3eV) 
because of a slight trigonal distortion of the material 
(see, e.g., LDA calculations with iVth order muffin-tin 
orbitals (NMTO) in Ref,^). As clearly stated in Ref,^, 
the interplay between strong electronic correlation and 



multi-orbital physics is expected to be the crucial in- 
gredient of the physics underlying the Mott MIT in in 
V2O3. In fact, the properties of the MIT in the Cr-doped 
V2O3 have been calculated (in some case even preceding 
the experimental results) by means of LDA-I-DMFT in 
Refs.— 1^, and, later, by including the orbital hybridiza- 
tion in Refi^. 

Beside the success in describing photoemission data, 
LDA-f DMFT can be also used to analyze optical spec- 
tra. While -at the DMFT level- the numerical effort for 
computing the optical conductivity 17(0;) is comparable to 
that for computing spectral functions, as vertex correc- 
tions can be usually neglectedi^, rough approximations 
have been always done in evaluating the optical dipole 
matrix elements in the localized (NMTO, Wannier, etc..) 
orbital basis. In particular, in the first LDA+DMFT cal- 
culations of (7 (a;) for V2O3S the dipole matrix elements 
were simply replaced by 1, while in later works^ the 
dipole matrix elements have been evaluated in the Peicrls 
approximation, including the effects of multiple atoms in 
the unit cell when necessarjiiSii^. 

Our results for V2O3 are summarized in Fig. |4l where 
we show in the first row LDA (left) and LDA-fDMFT 
(right) calculations for the optical conductivity obtained 
by using the optical matrix elements, while in the second 
row the corresponding calculations made with the PA are 
discussed (we use a = 4.95A , c = lAA as lattice param- 
eters, see Ref.— and references therein). In all cases, we 
also directly compare our theoretical results with the ex- 
perimental data reported in Refi^. Our analysis at the 
level of the optical conductivity clearly confirms the piv- 
otal role played by electronic correlations in the physics 
of V2O3: the LDA results show a much stronger Drude 
peak (almost an order of magnitude stronger) when com- 
pared to the experiments. The inclusion of correlations 
via DMFT significantly improves the situation: Due to 
the proximity to the Mott-Hubbard MIT one observes 
a marked spectral weight shifts from the Drude peak to 
higher frequencies, which makes the overall agreement 
with experiment much better in the region up to 1.5eV, 
where the experimental data are available. 

From our analysis, moreover, another important aspect 
emerges: in the case of V2O3 the PA (adopted in previous 
calculations, e.g.—) works satisfactorily well, at least in 
the low-energy t2g subspace: The improvements due to 
the inclusion of the full optical matrix elements only leads 
to small changes in the optical spectra up to 2cV both in 
the LDA and LDA+DMFT results, as it can be expected 
on the basis of the discussion of Section IIB, considering 
the small (or vanishing) value of the crystal field splitting 
between the localized t2n orbitals at the Fermi level. 



IV. CONCLUSION 

We have developed a program package for calculat- 
ing the optical conductivity on the basis of Wien2K, 
and make it available to the scientific community 



at www. wicn2k.at/reg-User/unsupportcd/wicn2wannicr. 
Electronic correlations, e.g., from DMFT, or finite life 
times, e.g., from impurity scattering, can be included 
via a corresponding self energy for the Wannier or- 
bitals. From this self-energy the Green function is cal- 
culated, which together with the full dipole matrix ele- 
ments yields the optical conductivity, disregarding vertex 
correctionsi^. 

The main topic of the paper is a careful comparison be- 
tween the dipole matrix element approach and the Peierls 
approximation, which is the de facto standard for lat- 
tice model calculations. We have considered two materi- 
als SrVOs and V2O3 as testbeds. The low frequency part 
(below 2-3 eV) of the optical conductivity stems from d-d 
transitions, at least for the two materials considered and 
many other transition metal oxides. This part is well cap- 
tured by the Peierls approximation. One can understand 
this by the high degree of localization of the degener- 
ate (or almost degenerate) Wannier d orbitals: Below ~ 
1 eV, for both vanadates, it is also sufficient to include 
only the three t2g bands out of the five d-orbitals. For 
the high frequency part (above 2 — 3 eV), on the other 
hand, not only the d Hubbard bands are relevant but also 
p-d transitions. This part of the spectrum is not well de- 
scribed by the Peierls approximation. Generalized Peierls 
approximation, while still approximate, also improves the 
description of p-d transitionsi^ at a computational cost 
comparable to the full dipole matrix calculation. 

The comparison to experiment shows that 
LDA-fDMFT with full dipole matrix elements gives a 
good description of the optical conductivity. In contrast, 
the Peierls approximation shows strong deviations at 
high frequencies. The same is true for the LDA optical 
conductivity even with the full dipole matrix elements. 
For instance, the LDA optical of SrV03 conductivity 
particularly shows a too pronounced peak at '^ 3.5eV. 
This peak stems from d-p transitions, and the DMFT 
correctly spreads the d orbital spectral over a larger 
energy region: Hubbard side bands are formed and the 
electron-electron scattering smears out the d bands. 

The residual differences between LDA-I-DMFT and ex- 
perimental infrared-spectra hence cannot be ascribed to 
the limitation of the Peierls approximation, but rather 
to effects beyond the LDA-I-DMFT scheme. For exam- 
ple, impurity scattering and the inclusion of non-local 
electronic correlations. The inclusion of the latter re- 
quires a considerable efforts of going beyond the stan- 
dard LDA+DMFT scheme, e.g., by cluster extension 
of DMFT—, dynamical vertex approximation^ or du- 
als fermion^, which also necessarily requires a proper 
treatment of vertex corrections. 
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